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Abstract 

We discuss the elementary physics of the final state Coulomb interactions in Hanbury- 
Brown Twiss interferometry, showing - with explicit comparison to E877 data for 7r + 7r~ 
and ir^p - that the Coulomb corrections in the pair correlation function can be well 
understood in terms of simple classical physics. We connect the classical picture with 
descriptions in terms of Coulomb wave functions, and investigate the influence of the 
"central" Coulomb potential on the pair correlation function. 

1 Introduction 

Hanbury-Brown Twiss interferometry of identical mesons has become an important probe 
of the evolving geometry of the collision volume in ultrarelativistic heavy-ion collisions [1-3] . 
The quantitative interpretation of the results depends critically on understanding the role 
of Coulomb interactions of the detected pairs of particles with each other, as well as the 
Coulomb interactions of the pair with the system of remaining particles. In this note we 
focus on the elementary physics of these processes, deferring detailed calculations to later 
publications. 

The simplest form of Coulomb correction is inclusion of the Gamow factor - the square 
of the relative Coulomb wave function, ^c(0), of the produced pair at zero separation - a 
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procedure which is followed in many analyses j|, H. The assumption is that the produced 
pair of identical particles is made in a relative Coulomb state (at zero separation), and 
the amplitude for doing so is thus modified from the bare amplitude by the factor ^c(O); 
non-relativistically 

where the dimensionless parameter n is given by 

zz'e 2 zz'a , , 

V= = r , 2 

Vrel V rel /C 

for a pair of particles of charges ze and z'e with relative momentum q = (p — p')/2 and 
relative velocity v re \ = q/m re d, with reduced mass m re d = m/2 for two particles of mass 
m. [We consider particles with \z\ — 1 here.] The Coulomb-corrected rate of production is 
then inferred to be the measured rate divided by |^<7(0)| 2 . For pairs of the same charge, the 
Gamow correction suppresses the probability of production at small q, by a factor tending 
at small q to 2Txr]e~ 2 ' Kri , while enhancing the probability of production for opposite sign pairs 
by a factor tending to 27r 1 77 1 at small q. 

Correcting for Coulomb effects by taking the relative Coulomb wave function at the 
origin is not physically correct in heavy ion collisions. As noted in Refs. || |7j (see also 
[8-11]), taking the finite size of the source size into account can produce significant effects. 
The standard Gamow correction assumes that the separation of the particles of the pair at 
creation is small compared with the (zero angular momentum) classical turning point, r t , 
defined by q 2 /2m re d = e 2 /r t . However, for pions, r t ~ (200 fm)/g 2 , where q is measured in 
MeV/c; for q ~ 10 MeV/c, a typical minimum value, r t is only 2 fm, and smaller for larger 
q. Since r t is much smaller than the characteristic heavy ion radius, most of the pairs of 
particles observed in a heavy ion collision are made at relative separations well outside their 
classical turning points. 

We note that for typical q, the three length scales in the problem - the turning point, 
the particle wavelength, and the two-particle Bohr radius, a = l/m re( ie 2 (= 387 fm for -kit 
and 222 fm for 7rp) - are cleanly separated: 

r t : 1/q : a = 2 : a q : (a q) 2 . (3) 

For 7T7T (or 7rp), a q = l/\r}\ = 1.96 (or 1.13) g/(Mev/c) ^> 1. The classical turning point 
is the relevant length scale here for Coulomb effects (not, as suggested in Refs. 0, [|, the 
two-particle Bohr radius). 

Furthermore, in the presence of many produced particles, the relative Coulomb interac- 
tion of a pair is highly screened, which also decreases effects of Coulomb suppression. The 
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motion of the particles in the pair is strongly affected by their interactions with the plasma 
of other particles, and the mutual Coulomb interaction of the pair becomes dominant only 
when the pair has sufficiently separated from the other particles in the system that there is 
small probability of finding other particles between the particles in the pair. (See Ref. |L2 



for a recent investigation of screening effects in the context of a particular model for high 
particle multiplicities.) 

The major effect of the Coulomb interaction between the particles of the pair, at distances 
large compared with r t , is to accelerate them relative to each other. Particles of the same 
charge are accelerated to larger relative momenta, thus depressing the observed distribution 
at small q, while particles of opposite charge are reduced in relative momentum in the final 
state, which builds up the distribution at small q. Although these effects are qualitatively 
similar those produced by the Gamow correction, they are quantitatively rather different. 

Our main focus in this paper is on a simple schematic model for effects of Coulomb inter- 
actions. As we shall see, correcting for the Coulomb final state interaction of non-identical 
pairs enables one to extract important information contained in the measured correlation 
function, about the spatial and temporal size of the emitting source. The same detailed in- 
formation about the Coulomb final state interaction is also important to correct correlation 
functions of identical particles for Coulomb effects. 

2 Toy model 

We construct a greatly simplified model to take the screening and acceleration effects into 
account by neglecting the Coulomb interaction between the pair for separations less than 
an initial radius r , and for separations greater than r including only the relative Coulomb 
interaction. Since the relative motion is in the classical region, conservation of energy of the 
pair implies that the final observed relative momentum q is related to the initial momentum 
of the pair q at r by (see, e.g., |5|, [13| 



_^ = ^L±£!, (4 ) 

2m red 2m red r 

where the upper sign is for particles of like charge, and the lower for particles of opposite 
charge. For example, for pions with r$ = 10 fm, q 2 = q" 2 ± 20(MeV/c) 2 . 

How does the acceleration affect the measured correlation function? In a heavy ion 
collision the distribution of singles is given in terms of particle creation and annihilation 
operators by 

n(p) = (at(p)a(p)), (5) 
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while the distribution of pairs of particles of momenta p and p' is given by 



nzijPiP') = ( a \p) a \p') a (p') a (p)) ■ (6) 
Interferometry experiments measure the pair correlation function 

cm = } n f:f}L m 

\n[p )n[p')\ 

where the braces in the numerator denote an average over the total momentum P for an 
ensemble of pairs from the same events at fixed relative momentum q = (p — p')/2, and 
in the denominator they denote an average over particle pairs drawn from different events. 
Note that the 4-vector product q ■ P vanishes identically, implying that the component of 
the relative momentum along the pair direction gy equals {E\ — E 2 )/2f3, where (3 is the pair 
velocity and (Ei — E 2 )/2 is the relative pair energy. 

Since the Coulomb interaction conserves particles and the total momentum of the pair, 
the final distribution n,2(p,p') of pairs of relative momenta q is thus given in terms of the 
initial distribution of pairs, n 2 (po ,Po'), by 

n 2 (p,p')d 3 q = n° 2 (p ,p ')d 3 q (8) 

Equation (f|), with changes in relative angles ignored, yields the familiar Jacobian (see, 
e.g., 0) d 3 q /d 3 q = qo/q. We can assume to good accuracy that the Coulomb interactions 
between pairs of particles negligibly affect the singles distributions in the denominator of 
(0). Then 

C(q) = q ^C (q ) = f 1 T 2 ^f] V2 C (q ). (9) 

To illustrate this toy model, we compare in Fig. 1 the predictions of Eq. @ with 
E877 data for the 7r + 7r - , n~p, and n + p systems produced in Au+Au collisions at the AGS 



14j, assuming that the bare correlation function Cq equals unity. Shown in this figure as 
dotted lines are the results of the toy model for r$ = 3 fm (rightmost curve), 9 fm, and 
15 fm (leftmost curve), along with standard Gamow correction (solid line). Except at very 
small relative momenta q £10 MeV/c, where effects due to the finite momentum resolution 
of the experiment become visible in the data, the model gives a good account of the data 
for tq in the range of 9 - 15 fm. By contrast, the Gamow factor considerably overpredicts 
the data for all q values shown here. As we see, the raw correlation data for non-identical 
particles contains information about the mean separation of pairs when screening effects 
become negligible, summarized in the toy model by the (possibly q dependent) parameter 
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r . Then we show the Coulomb correction factor deduced from Eq. ([5J) for tt + tc + in Fig. 2a 
and 7r~7r~ in Fig. 2b, for the same range of r$ (as in Fig. 1, the rightmost curve corresponds 
to ro = 3 fm). Again we see that use of the Gamow factor implies a correction which differs 
significantly from that of the toy model. 

With the initial radius r extracted from the unlike-sign data, one can then construct 
the Coulomb correction for like-sign particles^. Dividing the "raw" E877 data by the toy 
model correction factor, with r = 15/m, we obtain the correlation function for like-sign 
pions (crosses) shown in Fig. 3a for 7r + 7r + and Fig. 3b for tt~7t~, which also show the 
correlation function (vertical bars) derived by making the standard Gamow correction. Using 
the Gamow factor instead of the proper Coulomb correction leads to a correlation function 
which is ~ 30% wider, implying a correspondingly reduced radius parameter. Furthermore, 
the shape of the "Gamow-corrected" correlation function has considerable non-gaussian tails 
in the range 30 < q < 80 MeV/c. These tails do not exist in the raw correlation function 
and obscure the interpretation of the data. 

We note also that the procedure described here for Coulomb corrections is not restricted 
to one-dimensional correlation functions. Since the Coulomb correction depends only on 
the magnitude of the relative momentum of the pair and not on its orientation, one should, 
for example, in an mult i- dimensional analysis in terms of q out , q S ide, and qi ong , apply the 
correction for, say, each bin of q out , q S ide, Qiong separately, i.e., before projection onto the 
particular variable of interest. 

3 Connection with quantum-mechanical description 

The physics of the toy model above is contained in the Coulomb wave function describing 
the propagation of the pair. To recall how this works we suppress the total momentum of 
the pair, and focus only on the relative momentum. In the absence of Coulomb interactions 
the number of pairs of relative momentum q is given by 

N (q) = J d 3 rdVe^ ( ^)(J t (f)J(r')), (10) 

where J(r ) is the amplitude for creating a pair at separation r , and the brackets denote an 
average over the event. In the presence of Coulomb interactions we have instead 

N{q) = J d 3 rd 3 r'^c{rW c (r')(JHr)J(r')), (11) 

It is, in general, inadequate to correct the like-sign data with the inverse of the unlike-sign Coulomb 
correction. For 7r — n correlations this approximation is a rather good, except at very small q; it fails, 
however, for correlations among heavier particles. 
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where ipci?) = ipc{fyiFi{— Wi ^',i{qf — <f ■ r)) is the relative Coulomb wave function for a 
pair of relative momentum q at infinity. 

Pairs of low relative momentum have relatively low angular momentum, e.g., a pair 
produced at 10 fm separation with relative momentum 20 MeV/c can have at most one unit 
of relative angular momentum. Thus only the low partial wave components of the Coulomb 
wave function enter Eq. ([0]) with appreciable probability. We consider here just s-waves, 
for which we employ the WKB approximation to write the Coulomb wave function at radius 
r outside the classical turning point as0 

~ rp ( r )l/2 g l/2 Sill< ^ r )' ( 12 ) 

where the local relative momentum, measuring the rate of change of phase, 0, of the wave 
function, is given by 

1 10 

/ x d(p ( 2 2m red e 2 \ 

P{r) = -T- = \ q =F • 13 

dr \ r J 



(Equation (12), with ^-dependent (j)(r) holds as well for higher partial waves, £ > 0.) The 



normalization of (II2J) agrees with ([I]) as r — ► 0, while as r — > oo, the Coulomb wave function 
approaches 

t/j(r) — — sin(qr — rj In 2qr + 5q). (14) 
qr 

The results of the toy model follow if we assume that the source correlation function 
(J*(r)J(r')) is localized in both r and r' around r$. The correlation function for s-waves 
(denoted by superscript s), the absence of Coulomb, is 

N°(p) = [ d 3 rd 3 r' SinprSinpr ' (J^f)J(f% (15) 
J pr pr' 

then since in the region of any radius r outside the turning point the Coulomb wave function 
behaves locally as a free particle s-wave of momentum p(r), the correlation function is given 
by 

N s (q) « ^-N°(p(r )) (16) 

where the factor p(ro)/q arises from the denominators in Eq. (|1^) and (|15|). Consequently, 

C(q) ~ C (p(r ))p(r ) /q, the result in Eq. (|9]) with q = p(r ). 

2 The WKB approximation for the s-wave is quite good outside the collision volume for the parameters 
encountered in HBT interferometry in ultrarelativistic collisions. The condition for validity of the approxi- 
mation is \dp(r)/dr\ <c p(r) 2 , which for r <C clq, the region of interest, becomes the restriction, r^3/q 3 ^ 2 a^ 2 . 
For 7T7T (or 7rp) pairs with q > 20 MeV/c, WKB is reasonable for r down to ~ 5 fm (or ~ 6 fm). 
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With the connection between the toy model and the Coulomb wave function established 
we can now generalize the picture to extend to smaller values of the source radius Tq. In 
general, the effect of the Coulomb interactions depends on the detailed structure of the source 
correlation function (J^(f)J(f'))] we expect this correlation function to be approximately 
of the form 

(j\r)J(r')) « S{f , f')f{f' - f ); (17) 

S, which defines the spatial region of the source, varies as a function of its arguments on 
a scale of the size of the emitting region, and the Fourier transform of / defines the bare 
distribution of pairs, viz., 

« < W(£^''< ? >' (18) 

where, but for effects due to the finite size of the emitting region, N (q) ~ /((f)- 

For a first orientation we can identify the bare correlation of Eq. ( [TED with that exper- 
imentally observed for two-particle correlations in e + e~ annihilation [15-18], which indicate 
an HBT "radius" consistently below 1 fm. Thus we expect f(r — r ) to extend only over 
distances ^ 1 fm, or equivalently the bare correlation to vary on momentum scales of several 
hundred MeV/c. Since typical source sizes from an analysis of pion correlations following 
ultra-relativistic nucleus-nucleus collisions are of the order of 5 - 10 fm [|14], |3|, it is a reason- 
able approximation to neglect the size of the bare correlation and replace f{r — r ) in Eq. 
( |l7l) by a delta function. We write 

(f(f)J{f')) « 8(r - f')N S(r), (19) 

where S(r), of unit strength, describes the distribution of initial pair radii. With (|1^) we 
find 

C(q)/C (q) = N(q)/N = J dh\^ c {f)\ 2 S{r) (20) 

(an equation tracing back to Ref. ||). 

To illustrate the transition from the Gamow correction to the toy model we take S(r) 
to have a simple normalized gaussian form of range tq. S(r) = (27i)~ 3 / 2 rQ 3 exp(— r 2 /2rp). 
We show, in Fig. 4, for the 7r + 7r~ system, the results of calculations of C(q)/Co(q) using 



(|20|) for r = 1, 5, 9, and 18 fm (dash-dot curves, the highest for r = 1 fm, and falling 

with increasing r ). As r — * 0, the projection of the square of the Coulomb wave function 

onto the source S(r) converges to the standard Gamow correction (solid line).[] For larger r$ 

3 For tq < 0.1 fm (not shown in Fig. 4) the difference between the Gamow correction and a calculation 
with (|20|) is less than 0.5%. 
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values it rather quickly approaches the prediction of the toy model (shown here for an initial 
radius of 9 fm as a dotted curve), implying that, for pairs originating outside their classical 
turning point, the toy model provides an adequate and reasonably accurate description of 
the Coulomb effects. 



4 The "central" Coulomb potential 

We next turn to the question of the effects of the Coulomb interactions of the pair with the 
remaining particles. This is a difficult many-body problem, which we greatly simplify at 
this stage by assuming that the remaining particles can be described by a central Coulomb 
potential, Z e ffe 2 /r, where in a central collison of nucleus A with nucleus B the effective 
charge Z e fj is of order of the total initial nuclear charge (Za + Zb). This central potential 
accelerates positive mesons away and slows down the negatives, effects described by the 
Coulomb wave functions for the potential. The final momentum of any particle is related to 
the initial momentum p a at production point r a by 
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c(p) = c(p.) ± — • (21) 

where e(p) = (p 2 + m 2 ) 1 / 2 . (While Coulomb effects for the relative momentum can be 
treated non-relativistically as in Eq. (|J), the individual momenta are generally relativistic.) 
We ignore in this brief discussion quantum mechanical suppressions or enhancements of the 
amplitude for particle emission, as well as possible effects of angular changes in the individual 
particle orbits on the particle distributions. Then the single particle distribution is modified 
by the central potential, analogously to Eq. (§), by 

n(p) = n (p a )^y^ = Pae ^ n (p a ). (22) 
d A p p£{P) 

Both the magnitude of the distribution as well as its argument are shifted. 

Although the central potential shifts the singles distribution, it cannot introduce any 
correlations among emitted particles that have no initial correlation in the absence of the 
central potential, e.g., for different species or opposite charged pions, as one usually assumes. 
If in the absence of the central potential uncorrelated particles [C(q) = 1] are emitted in 
independent free particle states, then in the presence of the potential they are emitted in 
Coulomb states for the central potential, but still ri2(p,p') = n(p)n(p') and C(q) remains 
unity. 

For particles that are initially correlated as a consequence of Bose-Einstein statistics, 
n 2 (p,p') and n(p)n(p') will be modified both by the Jacobians of the transformations from 
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initial to final momenta, and shifts of argument. However, in forming C(q), the effects of 
the Jacobians in the numerator and denominator essentially cancel, and the primary effect 
is the shift in the arguments: 

C[q) {n(p a )n(p a >)Y (23) 

Since positive particles are accelerated, the final momentum difference, q = p — p ', of a 
positive pair will generally be larger in magnitude than it is initially, while for negative 
pairs the final momentum difference will generally be smaller. Thus we expect the central 
Coulomb potential to cause the size of collision volume extracted from positive pairs to be 
smaller than the actual size, and that from negative pairs larger than the actual size. As 
an illustration consider a pair of relativistic particles whose initial momenta p a and p a ' are 
equal in magnitude to p a , and final momenta p and p' equal in magnitude to p; then 

q = (P/Pa)q a ^ ?a (l ± ^"p/^ ) ■ (24) 

where the upper sign refers to both particles positively charged and the lower to both neg- 
atively charged. For Z ~ 150, r a ~ 7 fm and p a ~ 300 MeV/c, the effect is an increase 
for positives (and a decrease for negatives) in the observed scale of C(q) and decrease (or 
increase) in the extracted radius of ten percent. 

GB is grateful to S. Koonin for his classical wisdom, to J. Popp for helpful calculations, 
including those shown in Fig. 4, and T. Humanic and H. Strobele for discussions. This work 
was supported by NSF Grant PHY94-21309. 
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5 Figure captions 

FIG. 1. Comparison of the toy model, Eq. (||) for r$ = 3 fm (rightmost curve), 9 fm, and 
15 fm (leftmost curve), with E877 data [nj] for the systems 7r + 7r~, 7i~p, and vr + p, assuming 



a bare correlation function Cq — 1. Solid line: Gamow correction. 



FIG. 2. Coulomb correction, Eq. (^) for the systems tt + 7t + and tt 7r for the same range of 
To as in Fig. 1 (the rightmost curve corresponds to ro = 3 fm). 

FIG. 3. Toy model calculation of C(q) for like-sign pions (crosses), compared with the cor- 
relation function derived by making the standard Gamow correction (vertical bars). 
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FIG. 4. Transition from the toy model (dotted line, with r = 9 fm) to the Gamow correction 
(solid line) with decreasing source size, calculated from Eq. fl20|) (dash-dot curves). From 
highest to lowest dash-dot curves the source range ro is 1, 5, 9, 18 fm. The data shown are 

for 7T + /T~. 
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